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Detecting Sunyaev-ZeFdovich clusters with Planck: 
III. Properties of the expected SZ-cluster sample 
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1 INTRODUCTION 



The Sunyaev-Zel'dovich (SZ) effects jSunvaev & Zerdovichll972l 
1 198(1 1 Rephaeli 1995; Birkinshaw 1999) are promising tools for de- 
tecting clusters of galaxies out to high redshifts by their spectral 
imprint on the cosmic microwave background (CMB). This paper 
compiles the results of an extensive assessment of the observability 
of the SZ-effect for the European Planck-surveyor satellite based 
on numerical data, as described two preceeding papers (Schafer et 
al. 2004a,b). In the simulation, we try to model as many aspects of 
a survey of the CMB sky with Planck as possibly relevant to the 
search for Sunyaev-Zel'dovich clusters of galaxies and the extrac- 
tion of the weak SZ-signal. We discuss shortcomings of the sim- 
ulation and of the filtering schemes, quantify the properties of the 
resulting SZ-cluster sample and compare our results with previous 
studies. 

This paper contains a recapitulation of basic SZ-quantities 
(Sect. [2} and a brief outline of the simulation (Sect. |3j- The 
key results are compiled (Sect. |3} with a subsequent discussion 
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ABSTRACT 

The Planck mission is the most sensitive all-sky submillimetric mission currently being 
planned and prepared. Special emphasis is given to the observation of clusters of galaxies by 
their thermal Sunyaev-Zel'dovich (SZ) effect. In this work, the results of a simulation are pre- 
sented that combines all-sky maps of the thermal and kinetic SZ-effect with cosmic microwave 
background (CMB) fluctuations, Galactic foregrounds (synchrotron emission, thermal emis- 
sion from dust, free-free emission and rotational transitions of carbon monoxide molecules) 
and sub-millimetric emission from planets and asteroids of the Solar System. Observational 
issues, such as Planck's beam shapes, frequency response and spatially non-uniform instru- 
mental noise have been incorporated. Matched and scale-adaptive multi-frequency filtering 
schemes have been extended to spherical coordinates and are now applied to the data sets 
in order to isolate and amplify the weak thermal SZ-signal. The properties of the resulting 
SZ-cluster sample are characterised in detail: Apart from the number of clusters as a func- 
tion of cluster parameters such as redshift z and total mass M, the distribution n(<x)dcr of the 
detection significance <x, the number of detectable clusters in relation to the model cluster 
parameters entering the filter construction, the position accuracy of an SZ-detection and the 
cluster number density as a function of ecliptic latitude j3 is examined. 

Key words: galaxies: clusters: general, cosmology: cosmic microwave background, methods: 
numerical, space vehicles: Planck 



(Sect.|5J- The cosmological model assumed throughout is the stan- 
dard ACDM cosmology, with the parameter choices: fl M = 0.3, 
Q A = 0.7, H a = 100 ft km s~' Mpc~' with ft = 0.7, Cl B = 0.04, 
n s = 1 and erg = 0.9. 



2 THE SUNYAEV-ZEL'DOVICH EFFECTS 

The Sunyaev-Zel'dovich effects are very sensitive tools to observe 
clusters of galaxies out to large redshifts in submillimetric data. 
Inverse Compton scattering of CMB photons off electrons of the 
ionised ICM produces a modulation of the CMB spectrum and 
gives rise to surface brightness fluctuations of the CMB. One dis- 
tinguishes the thermal SZ-effect, where the thermal energy content 
of the ICM is tapped from the kinetic SZ-effect, where the CMB 
photons are coupled to the bulk motion of the ICM electrons, as 
opposed to the thermal motion of the electrons in the thermal SZ- 
effect. 

The relative change AT/T in thermodynamic CMB temper- 
ature at position as a function of dimensionless frequency x = 
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hv/(k B T C MB) due to the thermal SZ-effect is given by: 
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where the amplitude y of the thermal SZ-effect is the thermal 
Comptonisation parameter y. y is defined as the line-of-sight inte- 
gral of electron density times electron temperature. m c , c, k B and cr T 
denote electron mass, speed of light, Boltzmann's constant and the 
Thompson cross section, respectively. The kinetic SZ-effect arises 
due to the motion of the cluster parallel to the line-of-sight relative 
to the cosmological frame of reference given by the CMB: 
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IV is the radial component of the cluster's velocity. The amplitude 
w of the kinetic SZ-effect is called the kinetic Comptonisation. 

The SZ-observables are the line-of-sight Comptonisations in- 
tegrated over the cluster face. The quantities J/ and "W are conse- 
quently called the integrated thermal and kinetic Comptonisations, 
respectively: 
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(1a{z) denotes the angular diameter distance to the cluster at red- 
shift z. The fluxes of the thermal SZ-effect Sy(x) and of the kinetic 
SZ-effect S-w(x) as functions of observing frequency are given by 
eqns. ^6) and 0, respectively. The flux density of the CMB has a 
value of So = 22.9 Jy/arcmin 2 : 
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(6) 



(7) 



S^x) - ^ (expw _ i)2 

Exemplarily, table Q summarises the fluxes 5 j/ and Sw an d 
the corresponding changes in antenna temperature Ty and T-w for 
Comptonisations of J/ = "W = 1 arcmin 2 . 

The sensitivity of Planck is good enough to perform an SZ- 
survey of a significant part of the Hubble volume. Because the sim- 
ulation of such a vast volume including hydrodynamics resolving 
scales as small as cluster substructure is beyond the capabilities 
of current computers, a hybrid approach has been pursued for the 
map construction: All-sky maps of the SZ-sky were constructed 
by using the light-cone output of the Hubble-volume simulation 
IJenkins et all200lblColberg et all200(il) as a cluster catalogue and 
template clusters from the small-scale gas-dynamical simulation 
IWhite et alJ2002T) . In this way, sky-maps were constructed which 
contain all clusters above 5 x 10 13 Mq//i out to redshift z = 1.48. 
The maps show the right 2-point halo correlation function, incor- 
porate the evolution of the mass function and the correct distribu- 
tion of angular sizes. Furthermore, they exhibit cluster substructure 
and deviations from the ideal cluster scaling relations induced e.g. 
by the departure from spherical symmetry. The velocities used for 
computing the kinetic SZ-effect correspond to the ambient density 
field. The map construction process and the properties o f the re- 
sulting map are in detail described in Scha fer et alj 1200 4). Details 
of the thermal and kinetic SZ-maps can be seen in Figs. \T\ and [5] 
respectively. 
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Figure 1. Detail of the thermal Comptonisation map: A 2° X 2° wide cut- 
out centered on the ecliptic coordinates (/I,/?) = (100°, -50°) is shown. 
The smoothing imposed was a Gaussian kernel with A9 = 2'.Q (FWHM). 
The shading indicates the value of the thermal Comptonisation y, which 
is proportional to arsinh(10 6 y). The mesh size of the underlying Cartesian 
grid is ~ 14". 
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Figure 2. Detail of the kinetic Comptonisation map: A 2° X 2° wide cut- 
out is shown centered on the same position as Fig.0 i.e. at the ecliptic 
coordinates (A, ft) = (100°, -50°). The smoothing imposed was a Gaussian 
kernel with Ad = 2'.0 (FWHM). The kinetic Comptonisation w is indicated 
by the shading, being proportional to arsinh(10 6 w). 



3 SIMULATION OF SZ-OBSERVATIONS WITH PLANCK 

In this section, the simulation is outlined: First, the foreground 
emission components considered are summarised (Sect. |3~T1 . and 
instrumental issues connected to sub-millimetric observations with 
Planck are discussed (Sect. l3~2l . The data products resulting from 
the simulation at this point will be spherical harmonics expan- 
sion coefficients (S ( m ) v of the flux maps S v (6) for all nine ob- 
serving frequencies v, where the spectra have been convolved with 
Planck's frequency response windows and the spatial resolution of 
each channel is properly accounted for. Next, the signal extrac- 
tion methodology based on matched and scale-adaptive filtering 
is described (Sect. l3~3l . followed by the application to simulated 
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Planck channel 


1 


2 


3 


4 


5 


6 


7 


8 


9 


centre frequency vo 


30 GHz 


44 GHz 


70 GHz 


100 GHz 


143 GHz 


217 GHz 


353 GHz 


545 GHz 


857 GHz 


frequency window Av 


3.0 GHz 


4.4 GHz 


7.0 GHz 


16.7 GHz 


23.8 GHz 


36.2 GHz 


58.8 GHz 


90.7 GHz 


142.8 GHz 


resolution AO (FWHM) 


33'A 


26f 8 




9'.2 


7^1 


5f0 


5'.0 


5'.0 


5'.0 


noise level o"n 


1.01 mK 


0.49 mK 


0.29 mK 


5.67 mK 


4.89 mK 


6.05 mK 


6.80 mK 


3.08 mK 


4.49 mK 


thermal SZ-flux (Si/) 


-12.2 Jy 


-24.8 Jy 


-53.6 Jy 


-82.1 Jy 


-88.8 Jy 


-0.7 Jy 


146.0 Jy 


76.8 Jy 


5.4 Jy 


kinetic SZ-flux (S-w) 


6.2 Jy 


13.1 Jy 


30.6 Jy 


55.0 Jy 


86.9 Jy 


HO.OJy 


69.1 Jy 


15.0 Jy 


0.5 Jy 


antenna temperature ATy 


-440 nK 


-417 nK 


-356 nK 


-267 nK 


-141 nK 


-0.5 nK 


38 nK 


8.4 nK 


0.2 nK 


antenna temperature AT^y 


226 nK 


220 nK 


204 nK 


179 nK 


138 nK 


76 nK 


18 nK 


1.6 nK 


0.02 nK 



Table 1. Characteristics of Planck's LFZ-recefvers (column 1-3) and //F/-bolometers (column 4-9): centre frequency vo, frequency window Av as defined in 
eqns. |5J and 1101 . angular resolution A8 stated in FWHM, effective noise level <xn, fluxes (S y) and (Sw) generated by the respective Comptonisation of 
J/ = "W = 1 arcmin 2 and the corresponding changes in antenna temperature ATy and AT^y. 



Planck-data. (Sect. [3T4l . The morphology of peaks in the filtered 
maps as a function of signal profile model parameters is discussed 
(Sect. l3.5l and finally the algorithm for the extration of peaks in the 
filtered maps and the identification with objects in the cluster cat- 
alogue is described (Sect. l3~6l . A description of the software tools 
and the foreground modelling used for our simulation can be found 
in lReinecke et alj i2005t) . 

3.1 Foreground components 

We assume that the spectral properties of each emission component 
are isotropic, and that the amplitude of the emission is described 
by a suitably extrapolated template. It should be emphasised that 
this assumption is likely to be challenged, but the lack of observa- 
tion makes it difficult to employ a more realistic scheme. Further- 
more, it should be kept in mind that the extrapolation of the Galac- 
tic emission components by as much as three orders of magnitude 
in the case of synchrotron radiation is insecure. Details of the tem- 
plates and t he variou s emission laws are given in a precursor paper 
( Schit fer et all2004l) . 

• Cosmic microwave background: From the spectrum of Q- 
coefficients generated with CMBfast ISeliak & Zaldarriag3ll996h 
for the assumed ACDM-cosmology, a set of Of,,, -coefficients 
was derived by u sing the synalm code based on synfast by 
Iffivonet aljil99l . 

• Galactic synchrotron emission: The modelling of the galac- 
tic syn chrotron emission was b ased on an observation carried 
out by lHaslam et alj 1 198 lL 1 19821) at an observing frequency of 
v = 4 08 MHz, which has been adopted for Planck-usage by 
iGiardino et alj J2002I) . The spectral law for extrapolating the syn- 
chrotron flux incorporates a spectral break at v = 22 GHz, which 
has been reported by the WMAP satellite. 

• Galactic dust emission: The thermal emission from Galactic 
dust <SchlegeletaiTl 998) used a two temperature model proposed 
by C. Baccigalupi, which yields a good approximation to the model 
introduced by D.J. Schlegel. The emission is modelled by a super- 
position of two Planck-laws with temperatures T, = 9.4 K and 

= 16.2 K at a fixed ratio. 

• Galactic free-free emission: Modelling of the Galactic free- 
free emission was based on an //^-template provided by Finkbeiner 
(2003) and the spectral model proposed by IValls-Gabaud Jl998). 
which employs an approximate conversion from the H a intensity to 
the free-free intensity parameterised with the plasma temperature 
and describes the spectral dependence of the free-free brightness 
temperature with a v~ 2 -law. 



• Emission from carbon monoxide in giant molecular clouds: 
Rotational transitions of carbon monoxide give rise to a series of 
lines at integer multiples of the frequency Vc o _ 115 GHz. Its 
streng th is given by the template provided by iDame et"al] 1 19961 
2001). The line strengths of higher harmonics of the transition were 
determined by assuming thermal equilibrium of the molecular ro- 
tational states with an ambient temperature of T C o = 20 K. 

• Emission from planetary bodies of the solar system: We con- 
sidered a model of the heat budget of the planetary surface with 
periodic heat loading from the Sun and internal sources of thermal 
energy, heat emission and conduction to the planet's atmosphere 
and the back-scattering onto the surface by the atmosphere. From 
the orbital elements of the planet as well as from the motion of 
the Lagrange point L2 where Planck will be positioned we derived 
the distance to the planet and its position in order to compute a 
sky map. In total, we considered 1200 asteroids apart from the five 
(outer) planets, excluding Pluto which is to faint to be detected by 
Planck. 

The list of foreground emission components which possi- 
bly hamper with the observation of Sunyaev-Zel'dovich clusters 
is still not complete. Microwave point sources such as AGN and 
starforming galaxies were not included, likewise the modelling 
of zodiacal light originating from interplanetary dust in the Solar 
system was not attempted. Concise descriptions of Galactic fore- 
gro unds at C MB frequencie s derived from WMAP data are given 
by iBennett et all2003h and lPatanchon et alj <2004 . 

3.2 Instrumental issues 

The most important aspects related to the observations of the CMB 
sky by Planck are the properties of the optical system, the noise 
introduced by the receivers and the frequency response, which will 
be summarised in this section: 

• Planck-beam shapes: The beam shapes of Planck are approxi- 
mated by azimuthally symmetric Gaussians b(6) = exp ( — ^3 I 
with a-/) = The residuals from the idealised Gaussian shape 

Vo ln(2) 

are expected not to exceed the percent level. Table [J gives the an- 
gular resolution A9 for each channel. 

• Simulation of noise maps: Noise maps were generated by 
drawing Gaussian distributed random numbers from a distribution 
with zero mean and variance ctn given by TableQ These numbers 
correspond to the noise for a single observation of a pixel by a sin- 
gle detector. Consequently, this number is downweighted by y« dct 
(assuming Poissonian statistics), where n dcl denotes the number of 
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redundant receivers per channel, because they provide independent 
surveys of the microwave sky. In a second step, exposure maps were 
derived by simulating scan paths with Planck mission characteris- 
tics. Using the number of observations n obs per pixel, it is possible 
to scale down the noise amplitudes by yn„bs and to obtain a realis- 
tic noise map for each channel. 

• Frequency response and superposition of emission compo- 
nents: Adopting the approximation of isotropy of the emission 
component's spectral behaviour, the steps in constructing spherical 
harmonics expansion coefficients (S t m ) vo of the flux maps S(6, v) 
for all Planck channels consist of deriving the expansion coeffi- 
cients a tm of the template a(6), 



0(m 



dQ a{6) r;xoy «-» am = 2^2. afm Y "' m ' (8) 

f=0 m=-t 

converting the template amplitudes a (m to fluxes S( m , extrapo- 
late the fluxes with a known or assumed spectral emission law 
to Planck's observing frequencies, to finally convolve the spec- 
trum with Planck's, frequency response window for computing the 
spherical harmonics expansion coefficients of the average measured 
flux (S f m ) va at nominal frequency v by using: 



\S tm)y — 



J dv S rm (v)R V0 (v) 
fdvR V0 (v) 



(9) 



Here, S i m (y) describes the spectral dependence of the emis- 
sion component considered, and R Yo (v) the frequency response of 
Planck's receivers centered on the fiducial frequency vq. For all its 
channels, Planck's frequency response function R vo (v) is well ap- 
proximated by a top-hat function: 



Rvo(v) ■ 



1, v 6 [vo - Ay, Vq + Ay] 
0, v t [vo - Ay, vq + Ay] 



(10) 



The centre frequencies v and frequency windows Av for Planck's 
receivers are summarised in Table. Q In the final step, the averaged 
fluxes (S f,„) vo for each emission component are added to yield the 
expansion coefficients of the flux map. 

3.3 Construction of optimised filter kernels 

The basis of our signal extraction method is the concept of ma tched 
and scale-adaptive fi ltering pioneered by ISanz et all 1200 ll) and 
iHerranz et alJ 12002). which we generalised to the case of spher- 
ical data sets and spherical harmonics Yr„,(0) as the harmonic sys- 
tem replacing Fourier transforms and plane waves of the case of a 
flat geometry. The theory of matched and scale-adaptive filtering 
of multifrequency data sets is beautifully developed by the above 
mentioned authors in their formulation as the solution to a func- 
tional variation problem with increasing complexity. Filter kernels 
\jr(]0\ ,R) constructed for finding objects of a certain scale R modify 
the sky map w(0) by convolution and are required to minimise the 
variance cr 2 v (R) of the filtered map w(0, R) while fulfilling certain 
conditions: 

(i) There should exist a scale R for which the value of the fil- 
tered field (w(Ro)) at the position of a source is maximal, 

(ii) the filter should be unbiased, i.e. (w(Rq)) at the position of 
a source should be proportional to the amplitude of the underlying 
signal, and 

(iii) the variance cr 2 ,(R) should have a minimum at this scale Rq. 

In short, the filtered field should yield maximised signal-to-noise 
values for the peaks, while being linear in the signal and allow- 



ing the measurement of sizes of objects. We restrict ourselves to 
spherically symmetric source profiles superimposed on a fluctu- 
ating background which is a realisation of a homogeneous and 
isotropic Gaussian random field characterised by its power spec- 
trum. The filter resulting from the variation with the boundary con- 
ditions (i) and (ii) is called matched filter, and the filter obeying all 
three conditions is refered to as the scale-adaptive filter. 

Solving this variational problem in its extension to multifre- 
quency data and spherical maps yields spherical harmonics expan- 
sion coefficients i// r (£) for the filter kernels of each observational 
channel v as a function of the assumed profile of the source, of 
the spectral variation of the source flux with frequency and of the 
auto- and cross correlation power spectra C V(V .(/). The relative nor- 
malisation of the filter kernels in different channels encodes the 
optimised linear combination coefficients for co-adding the filtered 
maps. Apart from that, one obtains an analytical expression for the 
fluctuation amplitude a 2 = Y,fY,v i Yjv j >f / vj(£)Cv i vj{£)ilJvj{£) of the 
filtered and co-added maps, such that any peak height can be ex- 
pressed in units of the map's standard deviation, which directly al- 
lows to quantify whether a certain peak is likely to be a genuine 
signal or a mere background fluctuation. 

Filter kernels following from the matched and scale-adaptive 
multifrequency filtering algorithm were subjected to a thorough 
analysis. They are tested on two different data sets, one containing 
just CMB fluctuations and (non-isotropic) instrumental noise, and a 
second data set, which comprises all foregrounds in addition. From 
the comparison of the two data sets one will be able to quantify 
by how much the number of detections drop due to the foreground 
components and how uniform the cluster distribution will be pro- 
vided the removal of foregrounds can be done in efficiently. Details 
of the cross-channel correlation properties of the foregrounds as 
well as of the noise and the derivation of filter kernels and their 
properties are discussed in a precursing paper ISchafer et all2004l) . 



3.4 Filter construction and synthesis of likelihood maps 

Filter kernels optimised for detecting a family of King-profiles 
y{0) oc ^1 + (#/# c ) 2 ] were derived for a range of core radii 9 C and 
asymptotic slopes A. Specifically, seven values of 8 C , 

e c = o:o, 1:0, 2:0, 4:0, 8:0, 16:0, 32:0, (ii) 

and five values of A, 

A = 0.6, 0.8, 1.0, 1.2, 1.4 (12) 

were considered, keeping the large range in core radii in mind. 
Using different values for A is motivated by deviations from the 
generic baryonic profile and by asymmetric clusters. The sky maps 
were convolved with the filter kernels, co-added, normalised to unit 
variance, as described in the previous section and synthesised to 
yield likelihood maps. In the synthesis, all multipole coefficients 
up to £ = 4096 have been considered and the angular resolution of 
the resulting maps (A^idc = 1024, pixel side length ^ 3'4) is high 
enough to resolve single likelihood peaks. 

An important numerical issue of spherical harmonic trans- 
forms is the fact that the variance (measured in real space) of a map 
synthesised from the a [m -coefficients is systematically smaller with 
increasing C than the variance C(C) required by the a (m -coefficients 
on the scale A9 ^ n/£. This is compensated by an empirical func- 
tion, the so-called pixel window, which lifts the amplitudes a fm to- 
wards increasing values of I prior to the reconstruction. This effec- 
tively results in higher signal-to-noise ratios of the detected clus- 
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ters. In the numerical derivation of filter kernels very low multi- 
poles below t < 3 were excluded because of numerical instabilities 
in a matrix inversion and set artificially to zero. An important con- 
sequence of this will be discussed in Sect, l4.6l 



3.5 Morphology of SZ-clusters in filtered sky-maps 

Fig. [3] gives an impression how the morphology of a peak in the 
likelihood map changes if filter kernels optimised for the detec- 
tion of profiles with varying diameter and asymptotic behaviour 
are used. We picked an association of two clusters at a redshift of 
z — 0.1, which generates a signal strong enough to yield a signifi- 
cant detection irrespective of the choice of 6 C and A. 

The matched filter yields larger values for the detection signif- 
icance, which is defined to be the signal-to-noise ratio of the central 
object, in comparison to the scale-adaptive filter for that particular 
pair of clusters. Secondly, if filters optimised for large objects, i.e. 
large C and small A are used, the two peaks merge in the case of the 
matched filter, but stay separated in case of the scale-adaptive filter. 
Hence the scale-adaptive filter is more appropriate in the investiga- 
tion of neighbouring objects. Additionally, the matched filter seems 
to be more sensitive to the choice of 9 C and A. Within the range of 
these two parameters considered here, the significance of the clus- 
ter detection under consideration varies by a factor of four in the 
case of the matched filter, but changes only by 25% in the case of 
the scale-adaptive filter. 



3.6 Peak extraction and cluster identification 

It is an important point to notice that cluster positions derived from 
Planck are not very accurate. In this analysis, the SZ-clusters are 
extended themselves and possibly asymmetric, they are convolved 
with Planck's instrumental beams in the observation and recon- 
structed from filtered data, where an additional convolution with 
a kernel is carried out. Furthermore, the pixelisation is relatively 
coarse (typically a few arcmin). All these effects add up to a po- 
sition uncertainty of a few tens of arc minutes, depending on the 
filter kernel. 

All peaks above 3cr were extracted from the synthesised like- 
lihood maps and cross checked with a cluster catalogue. A peak 
was taken to be a detection of a cluster if its position did not de- 
viate more than 30!0 from the nominal cluster position. Peaks that 
did not have a counterpart with integrated Comptonisation i) larger 
than a predefined threshold value were registered as false detec- 
tions, likewise peaks were not considered that did not exceed the 
threshold value of 3<x in more than two contiguous pixels. In this 
way, a catalogue is obtained which is essentially free of false de- 
tections and where the fraction of unidentified peaks a mounts to 
5-7% for a realistic thr eshold of J/ min = 3x 10~ 4 arcmin 2 (Haehnelt 
1 1997tlBartelmanrl200 ll) . The cluster catalogues following from ob- 
servations with specific (8 C , /l)-pairs of parameters were merged to 
yield summary catalogues for both filter algorithms and both noise 
compositions. If more than one cluster is found in the aperture, the 
cluster with the largest value for the integrated Comptonisation is 
assumed to generate the signal. In the merging process, we deter- 
mine which choice of (8 C , A) yielded the most significant detection 
for a given object. 
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Figure 4. Distribution of pixel amplitudes d of the filtered and co-added 
maps, normalised to the variance <x predicted in the filter kernel deriva- 
tion, for a data set including CMB fluctuations and instrumental noise, fil- 
tered with the matched filter (upper left, solid line), for a data set including 
Galactic foregrounds in addition (upper right, dashed line), for a data set 
containing the CMB and instrumental noise, filtered with the scale-adaptive 
filter (lower left, dash-dotted line) and finally a data set with CMB, instru- 
mental noise and Galactic foregrounds, filtered with the scale-adaptive filter 
(lower right, dotted line). The filters have been optimised for the detection 
of beam-shaped profiles. 



4 RESULTS 

First of all, we investigate the noise properties of the likelihood 
maps (Sect, l4~TV followed by an analysis of the detection signif- 
icances (Sect. l4~2l the different filter algorithms are able to yield. 
Then, the number of detected clusters as a function of model pro- 
file parameters is investigated (Sect. l4~3l . The population of SZ- 
clusters in the mass-redshift plane (Sect. l4"4t . the distribution of 
the position accuracies (Sect. 14.51 and the spatial distribution of 
clusters (Sect. l4~6l are the main results of this work. Finally, the 
detectability of the kinematic SZ-effect is addressed (Sect, l4"7l , 



4.1 Noise in the filtered and co-added maps 

In this section, the statistical properties of the noise in the filtered 
maps is examined. The filter construction algorithm gives the vari- 
ance a of the filtered and co-added fields as a function of filter 
shapes iff vi {€) and cross-channel power spectra C VIY X€) by virtue of 
cr 1 = Yi(Hiljj'l'v i ('OCy. v .(i)i//y.(f)' Due to deviations from Gaus- 
sianity of many noise components considered (especially Galactic 
foregrounds), it is important to verify if the variance is still a sen- 
sible number. Fig. Ogives the distribution of pixel amplitudes for a 
combination of noise components and filtering schemes. 

Although the distribution of pixel amplitudes seems to fol- 
low a Gaussian distribution with zero mean and unit variance in 
all cases, there are notable deviations from this first impression. As 
summarised in Table|2| the mean of the distributions is compatible 
with zero in all cases, but the standard deviation is less than unity. 
Furthermore, the kurtosis of all distributions is nonzero, hence 
they are more outlier-prone as the normal distribution (barykurtic), 
which leads to a misestimation of statistical significances of peaks 
based on the assumption of unit variance of the filtered map, which 
the filtered map should have due to the renormalisation. This effect 
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Figure 3. An association of two clusters at z — 0.1, extracted with the matched multifilter (top row) and with the scale-adaptive multifilter (bottom row) from 
a map containing all Galactic compontents, CMB fluctuations and instrumental noise. The panel gives the likelihood maps and the statistical significances of 
the detection of the cluster at the image centre in units of cr, for varying values of 8 C with A fixed at A = 1.0. The side length of the panels is 4° . 



is strongest in the case of the matched filter. For the derivation of 
these numbers, only pixels with amplitudes smaller than \d\ < 4cr 
have been considered, such that the statistical quantities are domi- 
nated by the noise to be examined and not by the actual signal. The 
distributions are slightly skewed towards positive values, which is 
caused by weak signals below 4o\ The near-Gaussianity suggests 
that the residual noise in the filtered map is mostly caused by un- 
correlated pixel noise and filters seem to be well capable of sup- 
pressing unwanted foregrounds. 

Is it important to notice that the comparatively low threshold 
of 3cr imposed for extracting the peaks alone would yield a consid- 
erable number of false detections. This motivated the rather compli- 
cated algorithm outlined in Sect. 13.61 Supposing that the variance 
of the filtered maps is mainly caused by uncorrelated pixel noise 
which is smoothed to an angular scale of a 20' by the instrumental 
beam and by the filters causes the filtered maps to be composed of 
4n(180/n) 2 • 3 2 - 4 x 10 5 unconnected patches. Of these patches, a 
fraction of erfc(3/ v2) - 1(T 4 naturally fluctuates above the thresh- 
old of 3cr. In this way a total number of = 400 patches have signif- 
icances above 3cr. The requirement that the counterpart of the peak 
in the cluster catalogue generates a Comptonisation above a (con- 
servative) value of J/ m j„, i.e. that a cluster candidate is confirmed 
by spectroscopy, removes these false peaks from the data sample. 
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Figure 5. Distribution «(cr)dcr of the detection significances cr, for the 
matched filter (solid line, circles) in comparison to the scale-adaptive filter 
(dashed line, squares). The distributions are given for the clean data set in- 
cluding only the CMB, both SZ-effects and instrumental noise (thick lines, 
closed symbols) and in comparison, the data set where all Galactic fore- 
ground components are included in addition (thin lines, open symbols). 



4.2 Detection significances and total number of detections 

The distribution of detection sigificances is given in Fig.|5| One ob- 
tains about 10 3 detections at the significance threshold which drops 
to a few highly significant detections exceeding 20cr. At small cr, 
the scale-adaptive filter yields more detections than the matched 
filter, which catches up at roughtly 5cr. 

The total number of detections for each filter algorithm, for 
each data set and for two values of the minimally required Comp- 



tonisation J/ m i„ for spectroscopic confirmation are compiled in Ta- 
ble Due to its better yield of detections marginally above the 
threshold the scale-adaptive filter outperforms the matched filter 
by almost 30%. The reason for the increased number of low- 
significance detections is the systematically higher value of the 
variance of the residual noise field in the case of the scale-adaptive 
filter. The number of detections decreases by 25% if Galactic 
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filter algorithm data set mean n variance <x skewness j kurtosis k — 3 



matched 


CMB 


+ noise 




-0.0038 


+ 


0.0005 


0.9272 


+ 


0.0003 


0.0334 


0.5297 


matched 


CMB 


+ noise 


+ Galaxy 


-0.0009 


+ 


0.0005 


0.8902 


+ 


0.0003 


0.0154 


0.4232 


scale-adaptive 


CMB 


+ noise 




-0.0012 


± 


0.0005 


0.9090 


± 


0.0004 


0.0142 


0.2923 


scale-adaptive 


CMB 


+ noise 


+ Galaxy 


-0.0005 


± 


0.0005 


0.9023 


+ 


0.0004 


0.0076 


0.3125 



Table 2. Statistical properties of the filtered and co-added maps, derived from the first four moments of the amplitude distributions in Fig. [4] for all data sets 
and filter algorithms. The filters have been optimised for the detection of beam-shaped profiles. The errors given for the mean n and standard deviation cr of 
the distribution of pixel amplitudes correspond to 95% confidence intervals. 



foregrounds are included, relative to the data set containing only 
CMB fluctuations and instrumental noise. In a realistic observa- 
tion, one can expect a total number of ~ 6x 10 3 clusters of galaxies, 
compared to 8 x 10 3 clusters if only the CMB and instrumental 

noise were present. When comparing the total number of detections 

I 1 11 it I | III 
to analytic estimates (e.g. Aehanim et al. 1997; Kav et al. 2001; 

Bartelmann 2001), it is found that the number of clusters detected 
here is smaller, by a factor of less than two. 

One should keep in mind that the noise due to PLANCK's 
scanning paths is highly structured on the cluster scale and be- 
low, such that the assumption of isotropy of the noise is not valid. 
This has two important consequences: Firstly, assuming a simple 
flux threshold in analytic estimates is not valid because the noise 
is not uniform on the cluster scale and secondly the assumption 
of isotropy which is essential to the filter construction is violated 
which affects the sensitivity of the filters. 

4.3 Cluster detectability as a function of filter parameters 

The way the significance of a detection of a cluster changes when 
the core size 6 C and the asymptotic slope A are varied is illustrated 
in Fig.0for the matched filter. In general, the matched filter yields 
significances that are almost twice as large in comparison to the 
scale-adaptive filter for the specific example considered and con- 
sequently finds more clusters above a certain detection threshold. 
Furthermore, the matched filter shows a stronger dependence of 
the significance on the filter parameters C and A: The significance 
for the detection of the same object varies by a factor of four in 
case of the matched filter but only by 25% in the case of the scale- 
adaptive filter. This means that the derivation of cluster properties 
based on the filter parameter that yielded the most significant de- 
tection is likely to work for the matched filter, but not for the scale- 
adaptive filter. It should be emphasised, however, that the scale- 
adaptive filter keeps the likelihood distributions of the two objects 
from merging, in contrast to the matched filter. For that reason, the 
scale-adaptive filter may be better suited for the investigation of 
associations and pairs of SZ-clusters. 

Fig- El shows the number density of detectable clusters as a 
function of the King-profile's core size 8 C that entered the filter 
construction. Whereas the matched filter yields most detections at 
small values of 6 C , the scale-adaptive filter is better suited to de- 
tect extended objects. Most of the detections are registered at core 
sizes Q c = 8'. Additionally, the scale-adaptive filter's capability of 
detecting extended objects suffers from the inclusion of Galactic 
foregrounds, which cause the total number of detections to drop 
by 20%. In contrast, the matched filter is able to deliver a compa- 
rable performance for all values of Q c if Galactic foregrounds are 
included. 

The number density of clusters as a function of the King- 
profile's asymptotic slope A which the filters are optimised for is 
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Figure 6. Number density n(8 c ) of clusters as a function of the filter pa- 
rameter core size 6 C , for a data set including CMB fluctuations and instru- 
mental noise, filtered with the matched filter (circles, solid line), for a data 
set including Galactic foregrounds in addition (crosses, dashed line), for a 
data set containing the CMB and instrumental noise, filtered with the scale- 
adaptive filter (plus signs, dash-dotted line) and finally a data set with CMB, 
instrumental noise and Galactic foregrounds, filtered with the scale-adaptive 
filter (diamonds, dotted line). The thick and thin lines denote detections and 
peaks above 10~ 3 arcmin 2 and 3 X 10~ 4 arcmin 2 , respectively. 

given in Fig. The number of detections following from scale- 
adaptive filtering is relatively insensitive to particular choices of A, 
whereas the matched filter yields a higher number of detections in 
the case of compact objects, irrespective of the noise components 
included in the analysis. Fig.|8|illustrates how the number of detec- 
tions changes as a function of both C and A. It should be empha- 
sised that none of the graphs depicted in Figs. [6] 0and|8|shown the 
total number of detections for a given choice of 9 C and A and that 
graphs are not corrected for multiple detections of objects at more 
than (9 C , /!)-pair. 

4.4 Cluster population in the M-z-plane 

Scatter plots describing the population of detectable clusters in the 
mass-redshift-plane are shown in Fig [5] for the matched filter and 
in Fig. I10I for the scale-adaptive filter. The clusters populate the 
log(M)-z-plane in a fairly well defined region. There are only few 
detections beyond redshifts of z = 0.8, but the shape of the detec- 
tion criterion suggests the existence of a region of low-mass low- 
redshift clusters which should be detectable but which are not in- 
cluded in the map construction. It is difficult to predict the SZ prop- 
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filter algorithm 


data set 


J^miii = 10 3 arcmin 2 


y rain = 3 X 10 4 arcmin 2 


matched filter 


CMB + noise 


2402 


5376 


matched filter 


CMB + noise + Galaxy 


1801 


4199 


scale-adaptive filter 


CMB + noise 


3234 


8020 


scale-adaptive filter 


CMB + noise + Galaxy 


2428 


6270 



Table 3. Total number of detections in both data sets and with both filters, and for two values for minimally required Comptonisation. 




Figure 8. Number of detections n(6 c , A) as a function of both filter parameters core size 6 C and asymptotic slope A, for the matched filter (top row) in comparison 
to the scale-adaptive filter (bottom row). The figure compares the number density following from a clean data set containing the CMB, the SZ-effects and 
instrumental noise (left column) with a data set containing all Galactic components in addition (right column). n(8 c , A) is given for the minimal signal strength 
y min = 3 X 10~ 4 arcmin 2 (upper plane) compared to J/ m i n = 10~ 3 arcmin 2 (lower plane). 



erties of low-mass clusters because many complications in the sec- 
tor of baryonic physics come into play such as preheating, deviation 
from scaling laws and incomplete ionisation, which makes it diffi- 
cult to predict the number of clusters missing in our analysis. To- 
gether with K. Dolag we prepared an auxiliary SZ-map from a gas- 
dynamical constrained simulation of the local universe that would 
fill in the gap and provide clusters with masses M < 5 x 10 13 M o //i 
below redshifts of z < 0.1 jDolag et all2005l) . 

Fig. II ll gives the marginalised distribution in redshift z of the 
cluster sample. The shape of the redshift distribution is determined 
by the competition of two effects: With increasing redshift z the ob- 
served volume increases, but contrariwise, the number of massive 



clusters decreases as described by the Press-Schechter function and 
the SZ-signal becomes smaller proportional to d^ 2 (z)- Most of the 
clusters are observed at redshifts of z — 0.2 and the detection limit 
is reached at redshifts of z — 0.8. This applies to both filter algo- 
rithms and data sets alike. 

Fig. ^|gives the marginalised distribution of the cluster's log- 
arithmic mass m = \og(M/(MQ/h)). At high masses, both filtering 
schemes detect cluster reliably, but with decreasing mass, the filter 
algorithms start to show differences in their efficiency. The mass 
functions peak at a value of 2.5 x 10 14 M Q //i, and decrease towards 
smaller values for the mass due to the decrease in number density of 
objects and smaller SZ-signal strengh J/. Fig. ll3l gives the distribu- 
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Figure 7. Number density n(A) of clusters as a function of the filter pa- 
rameter asymptotic slope A, for a data set including CMB fluctuations and 
instrumental noise, filtered with the matched filter (circles, solid line), for 
a data set including Galactic foregrounds in addition (crosses, dashed line), 
for a data set containing the CMB and instrumental noise, filtered with the 
scale-adaptive filter (plus signs, dash-dotted line) and finally a data set with 
CMB, instrumental noise and Galactic foregrounds, filtered with the scale- 
adaptive filter (diamonds, dotted line). The thick and thin lines denote de- 
tections and peaks above 10~ 3 arcmin 2 and 3 X 10~ 4 arcmin 2 , respectively. 
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Figure 9. Population of clusters in the log(M)-z-plane detected with the 
matched multifilter for the data set containing the CMB, instrumental noise 
and all Galactic foregrounds. The minimal signal strength was required to 
be J/ m i n = 10~ 3 arcmin 2 . 

tion of the cluster's Compton-J/ parameter. The distribution is close 
to a p ower law as expected from virial estimates (c.f. ISchafer et alJ 
2004), but at low Comptonisations, all distributions evolve shal- 
lower, which is due to the fact that clusters fail to generate a peak 
in the likelihood map exceeding the threshold value. 

4.5 Position accuracies 

A histogram of the deviations between actual and reconstructed 
cluster position is given by Fig. 1141 The position accuracy is given 
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Figure 10. Population of clusters in the log(M)-z-plane detected with the 
scale-adaptive multifilter. Here, the detections are given for a data set con- 
taining the CMB, instrumental noise and all Galactic foregrounds. All peaks 
exceed aminimial Comptonisation ofJ/ m j n = 10~ 3 arcmin 2 . 
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Figure 11. Distribution n(z)dz of the detected clusters in redshift z, for the 
matched filter (solid line, circles) in comparison to the scale-adaptive filter 
(dashed line, squares). The figure compares detections in a clean data set 
containing the CMB, both SZ-effects and instrumental noise (thick lines, 
closed symbols) to a data set with all Galactic components in addition (thin 
lines, open symbols). Again iJ mm was set to 3 X 10~ 4 arcmin 2 . 



in terms of the squared angular distance A = 0^. because a uniform 
distribution would yield a flat histogram. The distribution is sharply 
peaked towards A = arcmin 2 . A fraction of 50% of all clusters 
are detected within 10' from the nominal source position, but there 
is a tail in the distribution towards larger angular separations. For 
most of the clusters, this position accuracy is good enough for di- 
rect follow-up studies at X-ray wavelengths, but not good enough 

for optical observa tions. 

As noticed in lHerranz et alJ 120021) . it is not trivial to assign 
a peak in the filtered co-added map to an actual underlying clus- 
ter. Apart from this complication, the simulation itself introduces 
many sources of positional uncertainties: The distance shown in 
Fig. I14l is the angular separation between the most-bound particle 
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Figure 12. Distribution n(m)dm of the detected clusters in logarithmic mass 
m = log(M/ (Mq IK)), for the matched filter (solid line, circles) in compari- 
son to the scale-adaptive filter (dashed line, squares). Here, the distributions 
are given for a data set including only the CMB, both SZ-effects and instru- 
mental noise (thick lines, closed symbols) in comparison to a data set con- 
taining moreover all Galactic foreground emission components (thin lines, 
open symbols). The minimal Comptonisation for spectroscopic confirma- 
tion was J/ m j„ = 3 X 10~ 4 arcmin 2 . 
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Figure 14. Distribution of the squared angular distance A = 6> 2 rc between 
actual and reconstructed source position on a great circle, for the matched 
filter (solid line, circles) in comparison to the scale-adaptive filter (dashed 
line, squares). The figure compares detections above 4cr (thin lines) with 
detections above 5<r for clusters detected with the parameter 8 C = S'.Q. The 
clusters were required to generate a Comptonisation J/ m j n exceeding 3 X 
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Figure 13. Distribution n(q)dq of the logarithmic integrated Comptonisa- 
tion, q = log(J/), for the matched filter (solid line, circles) in comparison 
to the scale-adaptive filter (dashed line, squares). Here, the distributions are 
given for a data set including only the CMB, both SZ-effects and instru- 
mental noise (thick lines, closed symbols) in comparison to a data set con- 
taining moreover all Galactic foreground emission components (thin lines, 
open symbols). 



in the Hubble-volume simulation and the peak in the filtered map, 
i.e. apart from the misalignment of the peak in the filtered map rel- 
ative to the maximum of the SZ-emission, there are misalignments 
of the SZ- template map relative to the cluster's barycenter and of 
the barycenter relative to the position of the most bound particle. 



4.6 Spatial distribution of Planck's SZ-cluster sample 

Fig. ll5l shows the number density of clusters as a function of eclip- 
tic latitude y = cos/3. The figure states that the Planck cluster 
sample extracted with the specific filters is highly non-uniform for 
low significance thresholds, where most of the clusters are detected 
on a belt around the celestial sphere, but gets increasingly more 
uniform with higher threshold values for the significance. This 
is due to the incomplete removal of low-£ modes in the filtered 
maps, which bea rs interesti ng analogi es to t he peak-background 
split iWhite et afll 19871 ICole & Kaiseul 19894) in biasing schemes 
for linking galaxy number densities to dark matter densities: Es- 
sentially, the likelihood maps are composed of a large number of 
small-scale fluctuations superimposed on a background exhibiting 
a large-scale modulation. In regions of increased amplitudes due to 
the long-wavelength mode one observes an enhanced abundance of 
peaks above a certain threshold and hence an enhanced abundance 
of detected objects. 

As Fig. ll6l indicates. the filtered and co-added maps do have 
large amplitudes for the hexadecupole which are certainly not in 
agreement with the near-Poissonian slope of C(l) oc I 1 typical for 
a random distribution of small sources. The incomplete removal 
of low-f modes shows that the assumptions about isotropy is vio- 
lated on large scales and C(€) ceases to be a fair description of the 
variance contained in the af,„-coefficients. Clearly, this is a serious 
limitation to the spherical harmonic approach. In general, the low-£ 
fluctuations are more pronounced for extended objects, i.e. large 8 C 
and small A, and they are stronger in the case of the matched filter 
compared to the scale-adaptive filter. 

Similarly, detection significances near the detection threshold 
are inaccurate due to the long- wavelength modes. A way to remedy 
this would be to introduce local estimates of the mean and variance, 
for example by considering the average and the standard deviation 
of the amplitudes in an aperture with a few degrees in radius. One 
must keep in mind that in the filtered map, the signal is strong and 
likely to affect these two values. It should be emphasised, however, 
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Figure 15. Number n(y)dy of clusters as a function of ecliptic latitude y = 
cos /?, for the matched filter (solid line) in comparison to the scale-adaptive 
filter (dashed line). The figure compares the number of detected clusters as 
a function of ecliptic latitude for detection significances > 4.2o" (circles, 
thin lines), > 4.8<x (squares, medium lines) and > 6. 0<t (diamonds, thick 
lines). 
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Figure 16. Power spectra C(C) of the filtered and co-added maps, where 
the filter kernels are derived for the parameters (9 C ,A) = (4'0, 1.0), for the 
matched filter and data set containing the CMB realisation and instrumental 
noise (red line), for the matched filter and the data set containing all fore- 
grounds in addition (green line), for the scale-adaptive filter and the data set 
including CMB fluctuations and instrumental noise (blue line) and for the 
scale-adaptive filter and the data set that contains all foregrounds in addition 
(yellow line). 

that the strongest signals exceeding values of = 6cr are uniformly 
distributed over the celestial sphere. 

4.7 Detectability of the kinetic SZ-effect 

In this section, we give the distribution of peculiar velocities in 
Planck's SZ-cluster sample, which is an important guide for ki- 
netic SZ-follow ups. As Fig.^]indicates, the distribution of pecu- 
liar velocities are well approximated by a Gaussian with zero mean 
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Figure 17. Number «(i/ p e C )d!/ pec of clusters, for the matched filter (solid 
line, circles) in comparison to the scale-adaptive filter (dashed line, 
squares). Again, the detections in a data set containing the CMB, both SZ- 
effects and instrumental noise (thick lines, closed symbols) are compared to 
a data set containing all Galactic foregrounds in addition (thin lines, open 
symbols). 

and standard deviation cr ve] = 300 km/s. For a dedicated search 
for the kinetic SZ-effect in Planck's SZ-cluster sample, velocities 
are drawn from this distribution, hence cluster bulk motions up to 
300 km/s can be expected in 68% of all cases and velocities in 
excess of 1000 km/s only for 11 to 16 objects, depending on the 
filtering scheme. 



5 SUMMARY AND DISCUSSION 

The properties of the likelihood maps and of the cluster catalogues 
following from applying matched and scale-adaptive filtering to 
the simulated flux maps are characterised in detail. According to 
our simulation, Planck can detect a number of = 6000 clusters of 
galaxies in a realistic observation with Galactic foregrounds (com- 
pared to over 8000 clusters if only the CMB and instrumental noise 
were present), which does not confirm the high numbers claimed 
by analytic estimates. 

• The noise properties of the filtered and co-added maps were 
examined in detail. It was found that the noise is very close to 
Gaussian after filtering, despite the fact that the initial flux maps 
had considerable anisotropic non-Gaussian features and despite the 
fact that the noise is highly structured and anisotropic on the clus- 
ter scale. Quantitatively, the variance of the filtered maps is smaller 
compared to the prediction based on the cross- and autocorrelation 
functions of the maps convolved with the filter. This discrepancy, 
which amounts to = 10% is due to numerics, but has the effect 
that significances of peaks are slightly underestimated.The clus- 
ter detectability as a function of filter parameters showed that the 
matched filter performs better on compact objects, where its deliv- 
ered significance depends strongly on the choice of A. The scale- 
adaptive filter works well on extended objects and is relatively in- 
sensitive to A. 

• The physical properties of the detected SZ-cluster sample 
made in terms of mass M, redshift z and integrated Comptonisa- 
tion J/: The cluster population in the mass-redshift plane is fairly 
well defined, and the marginalisation over the mass resulted in most 
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• Extragalactic point sources were excluded from the analysis 
due to poorly known spectra and clustering properties. In the sim- 
plest case of homogeneously distributed sources, there is a Poisson 
fluctuation in the number of point sources inside the beam area, 
which causes an additional noise component with power spectrum 
C(£) oc I 1 similar to uncorrelated pixel noise. If these sources have 
similar spectral properties, they could be efficiently suppressed by 
the linear combination of observations at different frequencies. 

• We did not attempt to simulate effects arising in the map mak- 
ing process and complications due to the 1 //-noise. So far it has not 
been investigated how well small structures can be reconstructed 
from time-ordered data streams. The map-making algorithms are 
chiefly optimised to yield good reconstructions of the CMB fluctua- 
tions by recursively minimising the noise, but to our knowledge the 
reconstruction of compact objects like SZ-clusters or minor planets 
has not been simulated for these algorithms. At the cluster scale, the 
dominating noise component is uncorrelated pixel noise, so that the 
contamination by 1 //-noise does not play a role on these scales. 

• Gaps in the data are a serious issue for the filtering schemes: 
Blank patches in the observed sky cause the power spectra C V , V2 (/) 
at different multipole order C to be coupled due to convolution with 
the sky window function. This is due to the fact that the Y fm {9, if>)- 
basis ceases to be an orthonormal system if the integration can not 
be carried out over the entire surface of the celestial sphere. Be- 
cause the linear combination coefficients are determined separately 
for each multipole moment i from the inverse of the covariance ma- 
trix CV, V2 (f), correlations between the covariance matrices at differ- 
ing t are likely to yield an insufficient reduction of foregrounds. 

• Galactic templates, especially the carbon monoxide map and 
the free-free map, are restricted to relatively low values in i and do 
not extend to high multipoles covered by Planck. For that reason, 
foreground subtraction at high values of i is likely to be more com- 
plicated in real data. Furthermore, one should keep in mind that the 
frequencies above 100 GHz are a yet uncharted territory and al- 
though the existence of an unknown Galactic emission component 
seems unlikely, the extrapolation of fluxes by two to three orders of 
magnitude in frequency may fail. 

The capability of Planck to detect SZ-cluste rs has been the 
subject of manv recent works, pursuing analvtical l Ashanim et al, 



of the clusters being detected at redshifts of z — 0.2, where the 
distribution starts decreasing to values of z — 0.8, where no clus- 
ters are detected. The distribution of detected SZ-clusters in mass 
M confirmed that the high-mass end of the Press-Schechter func- 
tion is well sampled, that most of the clusters detected have masses 
» 2.5 x 10 14 M o //? and that clusters of lower mass are increasingly 
difficult to detect. 

• The position accuracy is better than 10' in half of the cases, 
which is sufficient for X-ray follow-up studies, but the distribution 
exhibits a tail towards high discrepancies between the cluster posi- 
tion and the position of the peak in the likelihood map. 

• The investigation of the spatial distribution, especially in 
ecliptic latitude showed that the distribution of clusters gets increas- 
ingly uniform with increasing detection threshold. This is due to the 
fact that the filtered and co-added maps exhibit long-wavelength 
variations due to insufficient filtering at low multipoles. 

The simulation as presented here has a number of shortcom- 
ings that may affect SZ-predictions: 

• It was assumed for reasons of computational feasibiliy that all 
Galactic foregrounds had isotropic spectral properties. While this 
is an excellent approximation for the CMB, Galactic components 
can be expected to exhibit spatially varying spectral properties. For 
example, the spectral index of the Galactic synchrotron emission is 
likely to change with the propeties of the population of relativisitic 
electrons and the magnetic field and the spectrum of thermal dust 
changes with the dust temperature. The filter construction as it is 
would be applicable to those cases as well despite the fact that at 
fixed angular scale n/t, the cross power spectrum C VjV .(/) between 
frequencies v, and Vj ceases to be a good description of the cross- 
variance contained in the a;,„(v,)-coefficients. 

• Another related point worth mentioning is the approximate 
derivation of the covariance matrix from the a^v^-coefficients 
with the relation C YlV .(f) = (21 + 1)-' ZtL-e «<m(v,X m (v ; ) for the 
computation of matched and scale-adaptive filter kernels. This for- 
mula yields an unbiased estimate of the power spectrum C VjVj (i) 
only in the case of Gaussian random fields, which is certainly not 
the case if the Galaxy is included in the analysis. It would be more 
appropriate to derive the value of C v . r ,(f) that maximises the likeli- 
hood of describing the data, as described in.Bon d et al-Hl998l) . and 
to use this power spectrum for the derivation of filter kernels. 

• We did not include ICM physics beyond adiabaticity. Cooling 
processes in the centres of clusters give rise to cool cores, which 
can be shown to boost the line-of-sight Comptonisation y by a fac- 
tor of ~ 2-3. The volume fraction occupied by such a cool core is 
very small compared to the entire cluster and hence the total inte- 
grated Comptonisation J/ does not change significantly. For a low- 
resolution observatory like Planck, the primary observable is J/, 
and for that reason, SZ-observations carried out with Planck should 
not be affected by cool cores. A further complication is the ex- 
istence of non-thermal particle populations in the ICM, but their 
contribution to the SZ-flux modulation is very small. 

• There is a serious issue concerning completeness. The popula- 
tion of detections in the M-z plane suggests that low-mass clusters 
at redshifts z < 0. 1 should be detectable for Planck. This particular 
region of the M-z-plane is not covered by the SZ-map construction, 
but Planck would certainly add detections in this particular region 
of the Af-z-plane. The SZ-maps of the local universe provided by 
lHansen et alJ 120051) fill in this gap of low-mass, low-redshift sys- 
tems and can be combined with the SZ-maps covering the Hubble 
volume described in this work. 



| l997tlDelabrouille et alj2002tlBartelmanj200ll:lMoscardini et al 
2002) as well as semi - analytical lSjmz_gt_aU|200 0; Herran^yil 
2002t |Kav et alJl200lUDiego et alJl2002t iHobson & McLachlar 
2003; Geisbiisch et al. 2005) and numerical approaches l White 

• lAghanim et alJ J 19971) use analytical y6-profiles, an M-T rela- 
tion from /1-body data and the Press-Schechter function for gen- 
erating square SZ-maps with side length ~ 12°. To this map they 
superimposed CMB fluctuations, Galactic foregrounds and instru- 
mental noise. From this data they recovered the SZ-signal by multi- 
frequency Wiener-filtering outlined bv lBouchet et aljll999l) . They 
predict a total number of 7000 clusters with integrated Comptonisa- 
tions }/ t 9x 10~ 4 arcmin" and 10 4 objects at if i 5x 10~ 4 arcmin". 
These numbers sligh tly exceed our resu lts. 

• The paper bv iBartelmannl 1200 ll) and the related work by 
iMoscardini et alJ l2002t) take a purely analytical approach with 
spherically symmetric /J-profiles for describing the SZ-morphology 
and rely on the Press-Schechter function and an M-T relation from 
numerical data for predicting the SZ-signal of clusters. They in- 
corporate the effect of the finite instrumental resolution and re- 
quire the integrated Comptonisation if to exceed the value of 
3x 10~ 4 arcmin 2 . The total number of detectable clusters is stated to 
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be 10 4 , which again slightly exceeds our findings, but the distribu- 
tion of cluster masses M and the distribution of detectable clusters 
in redshift z is very similar to the results presented in this paper. 
The redshift distribution peaks at a very similar value, but extends 
to larger redshifts beyond z — 0.8. 

• T he papers written by ISanz et all 1200 ll) and iHerranz et alJ 
(2002), who developed the concept of matched and scale-adaptive 
multifiltering based on an extremal principle for flat topologies 
and Fourier-decomposition as the harmonic system, concentrate 
mainly on filter construction. They employ analytic SZ-profiles 
and describe the instrumental noise as uncorrelated Gaussian pixel 
noise, but c onsider the entire spectrum of Galactic foregrounds. 
iHerranz et ai]|2002h advocate a number of ^ 10 4 clusters to be 
detectable by Planck, which they estimate by extrapolating the av- 
erage number of detections in simulated 12° -wide patches to the en- 
tire celestial sphere, while restricting themselves to higher Galactic 
latitudes of \b\ > 19° (/ sky = 2/3). Similarlv. lpiego et aljEooll) ap- 
plies a Bayesian non-parametric method to the same data and finds 
a total number of 9 x 10 3 clusters at Galactic latitudes of \b\ t 12° 
(/sky = 0.8). Thus, both analyses are quite comparable with our 
approach concerning their estimated number of detections, while 
yielding cluster catalo gues tha t conta in slightly more entries. 

• In the study by Kav et al. (2001) the SZ-population was mod- 
elled using the Hubble-volume simulation as a cluster catalogue. 
The main aim is the difference of SZ-catalogues delivered by 
Planck in the ACDM cosmology compared to the tCDM model. 
The expected SZ-signal was derived based on an M-T-relation and 
they include an instrumental description including finite resolution 
and frequency response. By requiring a cluster to generate an in- 
tegrated Comptonisation exceeding the value of 3 x 10~ 4 arcmin 2 
in an area defined by the virial radius, they find a total number of 
5x 10 4 clusters of which a fraction of 1 percent is spatially resolved. 
The limiting redshift is stated to be z — 1.5, while the distribution 
in redshift peaks at a comparatively large value of z — 0.3 . . .0.4. 
In comparison, their cluster catalogue exceeds ours by a factor of a 
few, which is likely due to the fact that they neither consider instru- 
mental noise nor foregrounds, but concentrate rather on the fluctu- 
ating SZ-background alone as the mai n source of noise. 

• The Bayesian approach by Hobs on"& McLachlanl 1200 3l) fo- 
cuses mainly on the problem of peak finding and shows that 
the method they investigate can be readily applied to Planck- 
data. They consider a very simplified SZ-observation with Planck- 
characteristics at a single frequency, use analytic profiles, neglect 
all Galactic foregrounds and include only the fluctuating CMB and 
instrumental noise as noise sources. Consequently, they do not give 
astrophysical properties of the SZ-clusters their method is able to 
find, but it should be emphasised that the quantification of a peak 
height in terms of a Bayesian likelihood is far preferrable to our 
quantification in terms of a statistical significance. 

• In contrast , the filter scheme employed in the paper by 
iGeisbiisch et alJ 120051) is the powerful harmonic-sp ace maximum 
entropy method introduced by S tolvarov et aljj20 02). which is pri- 
marily optimised for component separation rather than the detec- 
tion of individual objects. The SZ-signal they put into the sim- 
ulation is determined from scaling relations and uses spherically 
symmetric analytic profiles. Including an accurate description of 
Planck's instrumentation and Galactic foregrounds, they find a total 
number of up to 1.1 x 10 4 ... 1.6x 10 4 clusters depending the M-T- 
relation for the choice cr 8 = 0.9, while these numbers decrease by 
^35% for a 20° -wide Galactic cut. Their distribution in redshift z is 
quite similiar in shape compared to ours - neither of us finds high- 
redshift clusters beyond z = 1 , but their distribution falls off slower 



with increasing redshift z. A grand result is their extraction of the 
SZ power spectrum, which our analysis due to its focus on the de- 
tection of individual peaks is not able to deliver. It should be kept 
in mind, however, that the component separation method, despite 
its prowess, assumes prior approximate knowledge of the emission 
component's power spectra, which are only partially available at 
HFI-frequencies above v = 100 GHz. 

• White 12003) puts emphasis on using a hydrodynamical sim- 
ulation of structure formation, albeit on a small scale, but includ- 
ing non-collapsed objects. He includes neither CMB fluctuations 
nor Galactic foregrounds, but chooses a high threshold value for a 
simple flux criterion for detection. The signal amplification strat- 
egy is smoothing with a Gaussian kernel and linear combination 
of maps. From the number of detections on a small patch of the 
sky he coarsely estimates the extrapolated number of 10 3 . . . 10 4 
clusters without giving details on their astrophysical characteris- 
tics, but concentrates rather on the observability in follow-up stud- 
ies. Extending this investigation, Schulz & White 1 2003) have con- 
sidered a range of future CMB experiments including Planck and 
modelled SZ-observations on the basis of a large rc-body simulation 
with a semianalytic receipe for the gas distribution inside dark mat- 
ter structures. They restrict themselves to high Galactic latitudes 
and include CMB fluctuations as well as (Gaussian) instrumental 
noise. Signal extraction is performed by applying a matched filter 
algorithm. They quantify catalogue completeness and emphasise 
complications in peak finding. 

In conclusion, the simulation presented in this paper demon- 
strates the abilities of Planck with respect to detecting Sunyaev- 
Zel'dovich clusters of galaxies even in the presence of anisotropic 
non-Gaussian noise components with complicated spectral depen- 
dences. Despite the fact that the high number of detections claimed 
by analytical estimates needs to be adjusted, it was shown that our 
results support the expectations on Planck's cluster sample and that 
the numerical tools for analysing the cross- and autocorrelation 
properties of all Planck channels and for filtering the data work 
reliably up to the high multipoles of I = 4096 considered here. The 
Planck catalogue of SZ-clusters of galaxies will s urpass X-ray cat- 
alogu es (e.g. the REFLEX catalogue compiled by Bohrinaer et al. 
120041 on the basis the Rosat all-sky survey) in numbers as it 
reaches deeper in redshift and is able to detect low-mass systems. 
It will contribute to the determination of cosmological parameters 
related to structure formation and dark energy, and shed light on 
baryonic physics inside clusters of galaxies. 
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